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We suggest two methods for the detection of the formation of disoriented chiral condensates in 
heavy ion colhsions. We show that the variance in the number of charged pious (in a suitable range 
of momentum space) provides a signature for the observation of a disoriented chiral condensate. 
The signal should be observable even if multiple domains of DxC form provided the average number 
of pious per domain is significantly larger than unity. The variance of the number charged pions 
alone provides a signal which can be used even if the number of neutral pions cannot be measured 
in a given detector. On the other hand, the probability distribution in R, the proportion of neutral 
pions to all pions emitted in heavy ion collisions in certain kinematic regions, has been suggested as 
', a signal of a disoriented chiral condensate. Here we note that the signature can be greatly enhanced 

■ by making suitable cuts in the data. In particular, we consider reducing the data set such that the k 

' pions with lowest pr are all neutral. We find that, given such cuts, (7?) can be substantially different 

from 1/3. For example, for a single T>xG domain without contamination due to incoherently emitted 
pions, (ii) is 3/5 given the pion with lowest pt is neutral, and 5/7 given the two pions with lowest 
Pt are both neutral, etc.. The effects of multi-domain D^C formation and noise due to incoherent 
. pion emission can be systematically incorporated. Potential applications to experiments and their 

' limitations are briefly discussed. 

m ■ 
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I. INTRODUCTION 

o 

00 . 

' During the past several years, there has been considerable excitement about the possibihty of the formation of 
^\ , disoriented chiral condensates (D^C) in heavy ion collisions [p]-p2|. The basic scenario is as follows: In an ultrarela- 
0^ ■ tivistic heavy ion collision some region thermalizes at a temperature above the chiral restoration temperature. If the 
system cools sufficiently rapidly back through the transition temperature, the region will remain in a chiral restored 
phase. However, this phase is unstable; small fluctuations in any chiral direction (cr, tt) will grow exponentially. This 
can create regions where the pion field has macroscopic occupation. It should be stressed that this scenario is not 
derivable directly from the underlying theory of QCD and contains a number of untested dynamical assumptions, 
principally that the cooling is rapid. Thus the failure of the system to form a DxC cannot be used to rule out that 
the system has reached the chiral restoration temperature. On the other hand, observation of the formation of a D^C 
would be clear evidence that the phase transition had been reached. 
K> , Unfortunately, since the scenario is not derived directly from a well-defined theory, it is difficult to know precisely 
5_j ■ what constitutes observation of a D^C. Assuming the system forms a single large domain of D^C containing a large 
[ number of pions there should be clear signatures. In the first place, one expects an excess in the number of low px 
pions produced. They would be at low pT since, by hypothesis, the region is large so the characteristic momentum 
is small; the excess would be measured relative to a purely statistical thermal distribution. Such a signal has the 
advantage of working even if multiple regions of D^C form provided each region is large enough so the characteristic 
momentum is sufficiently small to provide a discernible signal over the thermal background. Such a signal is not 
decisive since one could imagine some other collective low energy effects which produce low pT pions. 

A much stronger signal of a single large domain of D^C has been proposed. Since the pions formed in a D^C 
are essentially classical they form a coherent state. The coherent state has some orientation in isospin space (or 
more precisely the system is a quantum superposition of coherent states with different orientations and particular 
correlations to the isospin of the remainder of the system |l^,^l[]). In essence all of the pions in the domain are 
pointing in the same isospin direction. Provided the total number of pions in the domain is large, this implies that 
the distribution of the ratio of neutral to total pions in the domain is given by |^ P,pOt 

f{R) = (1.1) 

where R is the ratio of the number of tto's in the D^C divided by the total number of pions and /(i?) is the 



probability. The derivation of f{R) is quite simple and will be discussed below. The distribution in Eq. (1.1) is 
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qualitatively distinct from a purely statistical distribution in which the emission of charged and neutral pions is 
uncorrelated. The distribution from uncorrelated emissions in the infinite particle number limit approaches a delta 
function at R — 1/3. For finite (but large) particle number the statistical distribution is narrowly peaked about 1/3 
with a variance, (i?^) — (i?)^ = where Af is the total number of pions. Since these two distributions are so radically 
different one should in principle have a very clear signal if a single region of D^C where to form in heavy ion reactions 
and if the pions from the D^C are kinematically separated from other pions in the system. 

The dramatic nature of the preceding signature is based in large measure on the assumption that a single large 
domain of DxC is formed. A priori this seems rather unlikely for the following reason: If a large region of the 
system starts in a hot chirally restored phase and then rapidly cools through the phase transition, then there will be 
a large region which is unstable against growth of the pion field. Presumably, this happens as a "seed" fluctuation 
in a small region which rapidly grows. It takes a time of at least L/c for information about the formation of the 
domain to propagate a distance L. However during the time this fluctuation is growing out to L, the pion field at L 
has been sitting in an unstable situation. The characteristic time it can remain in this unstable configuration is r, 
the exponential growth time. If the information about the initial seed does not reach L is a time comparable to t 
the region near L will likely begin its own exponential growth but in a chiral direction uncorrelated from the initial 
growth. Thus, one expects domains of characteristic size cr Jl^ , p^ . 

The effect of multiple domains on the R distribution is fairly clear: it will tend to wash out the signal. If a 
large number of domains form and the pions emerging from different domains cannot be distinguished kinematically 
it is clear from the central limit theorem that the R distribution will approach a normal distribution. This normal 
distribution may be distinguished from the normal distribution arising from uncorrelated emission; the case of multiple 
domains of D^C will have a substantially larger variance. 

Unfortunately, there is an important practical limitation which makes it difficult to exploit the R distribution as 
a signature. Even under the most optimistic of scenarios, the total number of pions coming from D^C's will be a 
small fraction of the total number of pions. If one includes all pions produced in the reaction, the signal from the 
pions from the D^C will presumably be overwhelmed. Thus, it is highly desirable to use kinematic consideration to 
enhance the contributions coming from the D^C. In particular, it is sensible to study the R distribution for a sample 
restricted to low pT pions only. In any scenario where the D^C is well defined, i.e., the occupation number is large 
is likely to require a moderately large regions of D^C and the characteristic momentum spread in the D^C will be 
fixed by the inverse size of the region. Thus one expects D^Cs to preferentially produce moderately low px pions. 
(One also should restrict the pions in the distribution to a moderately narrow rapidity window). 

As an experimental matter, it should be relatively straightforward to cut on the momentum of the charged pions 
to select low pT pions in a given rapidity window. For neutral pions, however, it is not a simple matter. The neutral 
pions will decay in flight and will ultimately be detected as photons. If one is simply interested in the overall R 
distribution, without cuts, and if the detected photons come predominately from ir^ decays then one can use the 
n^/2 as a surrogate for n^o. Recent experimental searches have exploited this strategy p3|"p5[. However, in order 
to study the R distribution in a limited kinematical region it is necessary to reconstruct the tt* momenta from the 
observed photons in order to make kinematical cuts on the tt*' momenta. Since the number of neutral pions per 
event is large, the reconstruction of neutral pion momenta is likely to be a formidable task. This raises the following 
interesting question: Can one find a signature for tlie presence of regions of DxC of essentially the same quality as 
the R distribution but which does not require the measurement of neutral pions? 

On the other hand, even if one succeeds in reconstructing the neutral pions so that a low px cut can be applied, 
the noise may still be severe, as both the signal and the noise peak at i? = 1/3. So it is highly desirable to have some 
other cuts to further suppress the noise due to incoherently emitted pions. This raises another question: Can one 
construct additional cuts to further enhance the signal-to-noise ratio? 

In this paper, we will provide answers to both questions raised above. For the first question, we will show in Sec. II 
that the distribution of the number of charged pions (in a kinematically limited region) contains essentially the same 
information about DxC formation as the R distribution. This should greatly aid in searches for D^C formation. We 
also discuss additional information that can be inferred if n'^'s can be reconstructed. In particular, we show that 
the distribution of the total number of pions (in a limited kinematic region) provides a means to distinguish D^C 
formation from other hypothetical mechanisms for the production of a pion coherent state. Then in Sec. HI, we 
will suggest cuts which may dramatically enhance the signal-to-noise ratio. We study the conditional probability 
distribution of R given only for the events in which the k pions with the lowest px are all neutral, and we will show 
that the expectation value of R is shifted away from 1/3. Since incoherent emission will result in a very narrow peak 
around R = 1/3, any such shifts should be easily observable. Moreover, one can make successive cuts by increasing the 
value of k, and enhance the signal in each successive step. We believe these signatures will be useful in the searches 
for DxC formation in RHIC and LHC. 
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II. SIGNATURE OF DISORIENTED CHIRAL CONDENSATES FROM CHARGED PIONS 



In this section, we study the distribution of the number of charged pions produced from a D^C. We show that it 
contains essentially the same information about D^C formation as the distribution of R, and the distribution will be 
much wider in the presence of D^C formation than otherwise. This is of significance since it provides a signature 
which avoids the experimentally challenging task of reconstructing neutral pions. 



A. The Simple Case of a Single Domain 



We begin by studying an overly simple model and in subsequent subsections we will generalize our results to more 
realistic scenarios. In this simplified situation we assume that in every collision a single large domain of D^C is formed 
with a large particle number. Moreover, we will assume that the field strength and spatial distribution of this domain 
do not vary event by event, and that the pions produced in the D^C are kinematically completely distinguishable from 
all other pions in the system (including those pions produced from "a's" — i.e., fluctuations in the (qq) directions). 
Finally, we will assume that both isospin violating effects and explicit chiral symmetry breaking are negligible. 

By hypothesis, the region of D^C contains many particles and is essentially classical in nature. To simplify discussion 
we will adopt the usual convention of describing the physics in terms of the degrees of freedom in a linear sigma model 
with 0(4) symmetry, i.e., a and tt rather than directly in terms of the QCD degrees of freedom. We wish to stress, 
however, that we are not relying on the detailed dynamics of any particular variant of the a model. 

Since a D^C has macroscopic pion occupation numbers, it is tempting to describe it as a coherent state in the 0(4) 
chiral space as in Ref. However, in a generic heavy ion collision, the DxC (if produced) in the interior of the 

fireball will be quantum mechanically correlated with the high momentum emission at the edge of the fireball [p^ . 
As a result, the low px part of the physical state is not a pure state and is more appropriately described by a density 
matrix p. For example, when explicit chiral symmetry breaking due to non-zero quark masses are ignored, a chirally 
restored state would be described by the density matrix psym such that 

(a) - Tr(p,y„,a) = 0, (tt) - Tr(p,y„,7r) = 0; (2.1) 

and for the physical vacuum, 

(a) = Tr(pvac^T) - /, (^) - Tr(pvac^) = 0, (2.2) 

where the non-zero expectation value of a reflects that chiral symmetry is spontaneously broken. 

Now we are ready to study D^C states where the three angles describe the orientations of the D^Cs 

in the 0(4) chiral space. With the notation {X){t/j, 9, (p) = Tr(p£)^c'(?/;, 9, (p)X), where X is any operator, a D^C state 
is characterized by the following relations. 

{a)i^,e,(b)^Fcos^, 
{-kq){iP, 9, (p) — FsinV'COS^, 

(tTj.) (-(/;, 0, 0) = F sin V' sin cos (/), (2.3) 
{'Ky){tl}, 9,(j)) — F sin tp sin 9 sin (j), 

where = (cr)^ + (ttq)^ + {t^x)"^ + (ttj^)^ is independent of the angles -0, S and (j) but can depend on space and time. 
The number operators for neutral and charged pions are 

riQ ^ ttJtto, n± ^ tt^tTj, + tt^tt^, (2.4) 

and one can easily find their expectation values , i.e., {no (1^,0,4)) — TT{pDxc{i^,0,4')'>T'O,±)- It turns out that (no,±) 
can be factorized into the following form. 

{no,±){^,9,^)^{n)go,±{^,9,^), (2.5) 

where 

(n) - (no) (^ = 7r/2, = 0, <^) = (n±) (tA - 7r/2, = 7r/2, <^) (2.6) 

and (j) can take on arbitrary values. In other words, the expectation value (n) measures the total number of tt's 
produced by the D^C if fully oriented in a pionic direction. In general, (n) depends on the dynamical details of the 
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D^C, such as its probability distribution in position or momentum spaces. The geometrical factors go{ip,0,(l)) and 
5±('0j^7'/') will be called the neutral and charged proportions respectively, and they take the following forms. 

go{ip,e,(t>) ^ sin^ Tpcos'^ 9, g±{xl:,0,(j)) ^ sin^ ipshi^ 9. (2.7) 

Note that (n) does not depend on the orientation angles {i]j,9,<j}) while go_± do not depend on the dynamical details 
of the DxC. 



Equation (2.5) represents that both expectation values (no) and (n±) are related to (n), and hence to each other, 
up to geometrical factors go.±. Actually the relationship goes deeper than that: the probability distributions of the 
number of neutral and charged pion produced, denoted by PqIjiq) and p±{'n±) respectively, can both be expressed in 
terms of the same distribution p(ri). 

Po,±(?io,±;V', 6*, 0) =p("-o,±/ffo,±(V', 6*, 0))/3o,±(V', ^, 0), j dno^±po,±{no^±;ilj,9,(j)) ^ J dnp{n) = 1. (2.8) 
Notice that p{n) does not depend on the angles, and all the angula r de pendences in po,± come from the geometrical 



factors go,±- One can easily see that the above relation implies Eq. (2_^), and also the following relation on the higher 
moments: 

«±)(^,0,0) = (n'=)go^,±(^,0,0), (2.9) 

with 

{n'') = J n''dnp{n). (2.10) 

The probability distribution p{n) determines the distributions of neutral and charged pions produced from a D^C 
with a given orientation (as described by the angles {ip, 9, (j))) in the 0(4) chiral space. Let's recall that the probability 
distributions of the numbers of neutral or charged pion incoherently emitted in a heavy ion collision are Poisson- 
Gaussian, which have narrow peaks. To quantify the narrowness of a distribution, it is useful to define the deviance 
6[X] of a distribution of variable X such that 

(X^) = {l + S[X]){Xf, or {AX)^ = (X^) - (X)^ ^ S[X]{X)\ (2.11) 

For a Poisson-Gaussian distribution, 

(n^) = (n)^ + (n), (5[n] l/(n) ~> when (n) ^ oo. (2.12) 

If the distribution p{n) is wide (i.e., with non-zero deviance in the large (n) limit), the fluctuations in the numbers 
of coherently emitted neutral or charged pions will be much larger than the case of incoherent emission. As a result, 
substantial deviances in the number of neutral or charged pions produced may signify D^G formation. Unfortunately, 
the form of p{n) depends on the dynamical details of the D^G and is extremely model-dependent. Various forms 
of coherent states and density matrices have been motivated by different theoretical considerations [[l7| and lead to 



various forms of p(n). In particular, the well-studied Glauber coherent state |26 also give a Poisson-Gaussian p(n) 
and hence negligible deviances in the pion numbers in the limit (n) — > oo. So it may seem impossible to determine 
whether D^Gs are formed by studying the fluctuations in the number of neutral or charged pions produced. However, 
recall that the above discussion presumes we have a D^C with a definite orientation in chiral space, i.e., definite 
values of tp, 9 and (j)- In DxO formation, however, the orientation is randomly generated, and the relevant object of 
study is an density matrix with some appropriate probability distribution functions for the angles, which we will now 
proceed to study. 

The above analysis shows that, for each set of orientation angles {'ip,9,(j)), the distributions of no,± are normal. 
However, since the orientation is randomly generated in the process of spontaneous symmetry breaking, one cannot 
predict {ip, 9, (p). On the other hand, since we are neglecting explicit chiral symmetry breaking, the system is equally 
probable to point in any direction in chiral space. Moreover, since by hypothesis we are in a semiclassical situation 
(large (n)), it is legitimate to work with probabilities rather than amplitudes. Using the technology of Ref. |2^, the 
probability distribution in the angular variables is given by the unit measure: 

d^P{^,9,(t)) = -i^ sin^i/; sin 9 di/j d9 dcj), (2.13) 
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and the relevant density matrix is 



PDxC 



(2.14) 



One can use Eq. (2.7) to reparametrize the probabihty distribution (2.13) in terms of the neutral and charged 
proportions go.i — {^o,±) /{n): the fraction of neutral or charged pions among all particles produced by the D^C. 



(2.15) 



^ V5o(l - go - 9± 

From this one can obtain the marginal probability distribution of and g± by integrating over the other variable. 



dP[go) = foigo)dgo = - \ - 

TT V 50 

dP{g±) = f±{g±)dg± = dg±. 



go 



dgo, 



(2.16) 



These distribution functions are plotted in Fig. 1 in solid curves. It is obvious the both distributions are far from being 
normal. The function /o is heavily skewed towards the low end and actually diverges as 1/ y/go when go — > 0. On the 
other hand, f± is flat, and g± is equally likely to assume any value between and 1. This is drastically different from 
pion emission from an uncorrelated source (the dotted curves in Fig. 1 1^), where both distributions would be normal. 




More quantitatively, one can calculate the first and second moments of go and g±. 

{go) = 1/4, (<?2) = 1/8, 

{g±) - 1/2, (gl) = 1/3, (2.17) 

and hence the respective deviance 5[go] and (5[g±], 

S[go] - 1, S[g±] = 1/3. (2.18) 

What does the distributions of the proportions go,± tell us about the distributions of 7io,±? It is obvious that when 
n is fixed, (?o,± give the pion distribution. In reality, of course, n is not fixed; we have shown that it behaves like a 



^ For uncorrelated emissions, the probability distributions are Poisson-Gaussian with mean 1/4 and 1/2 for neutral and 
charged pions respectively. The variances depend on the number of independently emitted pions; the plots correspond to the 
case of n = 50. 
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Poisson distribution if we have a conventional Glauber coherent state p^ . However, since Poisson distributions are 
sharply peaked, the dispersion of n will simply smear the distribution of ^^o.tt 

slightly, without changing the overall 

shape qualitatively. More specifically, note that |^ 

(no,±) = (n) (.go,±), 

K±) - in') (5o' ±> = (1 + '^N) (1 + S[go,±]) {{n) (5o,±))' ^ (1 + '5K±])K±)^ (2.19) 



The first equality gives 



(no) = (n)/4, (n±) ^ (n)/2, (2.20) 



which have the simple interpretation that, given the symmetry between the four directions (tTx, ttq, a) in the chiral 
space, a quarter of the particles produced by the D^C will be ttq, while half of them will be tt^ or TTy. On the other 
hand, the second equality gives 

S[no,±] = S[n] + S[go^±] + S[n]S[go,±] > S[go,±]. (2.21) 

For a Glauber state, 5[n] = l/{n), S[go] = 1 and (5[g±], one has 

5[no] = l + 2/{n), (5[n±] = 1/3 + 4/3(n). (2.22) 

As expected, when (n) is large, the deviances i5[no.±] approach (5[(7o,±]- For other states with non-zero S[n\, the 
deviances <5[no,±] would be even larger. Note that both (5[7io,±] are of order 1, in contrast to an uncorrelated emission 
which would have S = \/ {n). In an ideal world, such enhancements of J's would indicate the existence of D^C. 
It is interesting to compare our analysis with that in Ref. [l^ , where (5[no,±] (called Cm.+ -+- in the reference) has 



also been calculated with (5[no] = 4/5 and 5[n±] = 1/5, in contrast to our results in Eq. ( 2.18 ) with the values 1 and 
1/3 respectively. The discrepancy originates from a key difference in the physical scenarios: the states in Ref. |l^ is 
equally likely to point at any direction in the 0(3) isospace (without a a direction) while our state is equally likely to 
point at any direction in the 0(4) chiral space. The extra uncertainty of the orientation of the D^C in the a direction 
is reflected in the extra angle ip in addition to the usual polar and azimuthal angles 6 and (j) of the sphere in isospace. 
Even with a fixed direction in isospace and hence fixed 9 and varying ip from to 7r/2 will change the pion number 
of both species from zero to its maximum value. This extra uncertainty in the a direction increases the fluctuation 
in the number of pion produced, leading to larger (5[no,±] in our analysis than their counterparts in Ref. [ p^ . By 
hypothesis, all direction in the 0(4) chiral space with the a direction are equally likely, and we believe our formalism 
is more appropriate than that in Ref. p2[ in describing a D^C. 



B. More Realistic Scenarios: Multidomain Formation 



We have studied the simple case of pion emission from a single large (in the sense that (n) is large) D^C domain. 
As discussed in the introduction, this scenario is presumably not realistic. We will now proceed to study more realistic 
scenarios with multi-domain formation. The main point is, even though the probability distribution is smeared out 
because of the lack of alignment (in the chiral space) between the different domains, one feature survives, namely 
the large variance of the distributions. In particular, we will see that the variances for both neutral and charged 
productions are still much larger than that of an uncorrelated emission. 

Let's consider the case which we have domains, each with the same (n). We are also assuming both (n) and N 
are much larger than unity. The total number of pions of each species produced is the sum of pions of that particular 
species produced in each domain, the distribution of which has been discussed in the previous section. 

N N 

&10 = '^o ^ ' E "± ■ (2.23) 

j=l i=l 



^ The two set of brackets in the right-hand side of the equations below have different physical origins. The expectations 
of (7o,± are statistical in nature, while that of n is via quantum mechanical smearing of a coherent state. Consequently the 
distributions of n and go,± are assumed to be uncorrelated, and we derived the equalities below. 
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By the central limit theorem, the probability distribution of SriQ.i will approach normal distributions when N is 
large. However, we will see that the variances of the Gaussian distributions will be much larger for pion production 
from a D^C than those of uncorrelated pion emission. 

Since the pion production in each domain are independent, jt-q j- are independent random variables. Hence the 
mean of Sno,± is just the sum of the means of all n-g'j-i 

N 

(Eno,±) = - iV(no^±) = Uo,±, (2.24) 

2 = 1 

(i) 

and the variance of the sum no,± is just the sum of the variances of each 

N 
i=l 

= N{no,±)H[no,±] = {S[no,±]/N)Ml±. (2.25) 

In other words, 

Sno,± — = — — = . (2.26) 

In comparison with uncorrelated pion production, with 6 = l/A/'o,± we see that for multi-domain D^C the deviances 
are enhanced by a factor of eo,± = <5[^o,±](?^o,±)- 

£0 = ((n) + 2)/4, e± = ((n) + 4)/6. (2.27) 

A priori (n) can take any value, but eo,± are larger than unity for any value of (n) > 2. Even for a very modest 
(n) = 8, Co = 2.5 and e± = 2, leading to substantial widening of the corresponding distributions, an observable 
signature of D^C formation. For larger values of (n), the broadening will be even more pronounced. 

While the above scenario describes D^C with multi-domains, a probable feature of D^C formation in the real world 
(if it happens at all), it is still unrealistic in assuming all the domains are of equal strength, i.e., with the same (n). 
Instead one expects (n) of different domains to fall under a certain probability distribution, which depends on the 
details of the model. Naturally, one questions if the signatures discussed above s till su rvive under such circumstances. 

Let's consider the case with N domains, with different (n^*^) ^ 1. Equation ( ^.24| ) becomes 

N 

(Sno,±) - = ^Kiy = A/'o.±, (2.28) 

1=1 

where (rio.±) is the average of (ri,Q*j_). Equation (2.25) becomes 

N N 

(ASno,±)^ = E(^"o:L)' = E(-o:L>'5[n«±]- (2-29) 

i=l 1=1 

By the inequalities S[n^^U > '^[5o.±] Eq-dUl)) and (">^ > (Ef (">)V^ (mean of squares is larger than 
square of mean), we have 

(ASno,±)2 > (E(4:±)') S[go,±] > (E(4:±>) S[go.,±]/N ^ iS[go.,±]/N)Ul^. (2.30) 
In other words, 



cr^ 1 ASno,±)^ S[gQ^±] S[go^±]{no^±) 

d Sno,± = — > — — = . (2.31) 
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(Compare Eq. (2.26).) Again, the deviances are much larger than that of uncorrelated emission with 5 = l/A/o,± 
when (n) ^ 1 by the foUowing enhancement factors. 



£0 > (ri)/4, e± > {n)/6. (2.32) 

So we see that, even with domains of unequal strengths, the number of neutral or charged pions produced by a D^C 
will still have a much wider distribution than that from independent, uncorrelated emission. 

Lastly, one may also ask if the distribution of Sno,± will approach normal distributions when N is large in the 
case of domains with unequal strengths. In this case n^*-' 's do not all fall under the same probability distribution and 
the most simple form of the central limit theorem does not apply. On the other hand, there are generalized forms 
of the central limit theorem, which state that as long as the probability distributions are sufficiently "well behaved" , 
the sum of N random variables will still fall under a normal distribution when iV — > oo. It is actually possible to 
argue that the distribution of Sn± does approach a normal distribution by the Lindeberg generalization of the central 
limit theorem. (See, for example. Sec. 6.E of Ref. [^.) Whether the same conclusion holds for Eng is still an open 
question. 

Let us recapitulate what we have shown: one can calculate the probability distribution of the number of neutral or 
charged pions produced as a result of DxC formation. The resultant deviances (5's are not of the order 1/Af as in an 
uncorrelated emission, but are instead enhanced by factors e's which are of order (n). Seeing such enhancements of 
deviances would be signatures of coherent pion productions. 

One can understand the origin of such enhancements of statistical fluctuations of the number of neutral or charged 
pions by considering the following analogy. Consider two groups of gamblers playing roulette in a casino: N lawyers 
at $100 tables, and lOOA^ physicists at the $1 tables, where the odds are the same. If each lawyer and physicist is 
given n chips, of $100 and $1, respectively (so that the total amount given to the lawyers, J\f = N x lOOn, is the 
same as that given to the physicists, JV = lOOA^ x n), and is required to bet all of them, the average loss will be the 
same for both groups as long as they are following the same betting strategies. However, it is easy to see that the 
statistical fluctuation of the loss of the lawyers would be much larger for that of the physicists. In other words, the 
amount of loss, as well as its standard deviation, is "quantized" in units of the value of the bets. The larger the bet, 
the larger the fluctuation. On the other hand, one can also turn the argument around; a discerning external observer 
can deduce, with the knowledge of the gambling strategies, the size of the bets of the lawyers from the statistical 
fluctuations of the lawyers' losses, and do likewise for the physicists as well. 

Just as the chips in a single bet of a lawyer share the same fate (either win or lose), all the pions in a single D^C 
domain share the same orientation in the chiral space. As a result, the fluctuation of the number of pions in each 
direction in the chiral space is enhanced by a factor proportonal to (n), the number of pions in each D^C domain. 
And by reverse argument, one can deduce whether coherent pion emission is taking place by measuring the fluctuation 
of the number of emitted pions. 



C. Including Incoherent Contributions 



Our analysis has assumed that all the pions originate from coherent emissions, and each of these coherent states 
has (n^')) '3> 1. In the real world, there is contamination from independent pion emissions, and the total numbers of 
neutral or charged pions are the sums of these two contributions. 

no,±-41+4"i^\ (2.33) 

where "c" and "inc" stands for "coherent" and "incoherent" , respectively. The variances of the numbers of coherently 
produced neutral or charged pions are enhanced while those of incoherent production arc not. 

(A4:i)^ ^ £o,±(4:L>, (AnS±^)^ = (<±^>, (2.34) 



where eo.± are of order (71) ^ 1 (cf. Eq. ( 2.27 ), ( 2.32 )). Then one can calculate the variance of the sum of the two 
contributions. 

{Anf = {A4lr + (An(;"±^V = 2o,±(«), ~^o.± = X6 + (1 - x) = 1 + x(e - 1), (2.35) 

with 

X = {ni%)/M, l-x= V("o,±). (2.36) 
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The parameter x measures the fraction of pions which are coherently produced: x = 1 when all the pions are from 
D^C, while x = when all of them are independently emitted. Obviously, the more incoherent pions in the sample, 
the smaller is the enhancement factor eo,±- 

While these incoherently emitted pions dilute our signatures for D^C, they have different momentum spectra from 
those from DxC. DxC pions, being produced from coherent state, carry low pT- The typical is of the order of l/L, 
where L is the size of the domain from which the pion originates. In contrast incoherently emitted pions can carry 
high pt- Therefore applying a low px cut can minimize the noise from incoherent pion emissions. 

It is also advantageous to measure the rapidity of the pions and count their numbers in narrow rapidity windows. 
Bear in mind that the rapidities of the pions are, up to small dispersions, equal to that of the original domain. As we 
have mentioned, it is probable that many domains are formed in a single collision, and all these domains may have 
different rapidities. For example, the domains at the surface of the "fire ball" are moving with high speed relative to 
the domains at the center. By binning the pions according to their rapidities, one can partially separate the pions 
from different domains, and the signals are enhanced as a result. 



In summary, we suggest the following procedure in looking for signatures of DxCs. 

• Count the number of neutral or charged pions event by event from heavy ion collision experiments and measure 
their individual transverse momenta and rapidities. 

• Apply a low pr cut to suppress the noise due to uncorrelated pion emission. 

• Bin the events in different rapidity windows. 

• In each rapidity window, plot the number of events vs. the number of neutral or charged pions in histograms. 

• Evaluate the mean, (rio,±), and the variance, (Ano,±)^, in each rapidity window. 

• If we find (AriQ.i)^ is substantially larger than (no_±), then we are seeing possible signatures from DxCs. 

The above procedure allows us to search for signatures from DxCs by counting only the charged pions. This is 
important as, with our present technology, it is difficult to count the number of tto's in a momentum bin, which would 
mean reconstructing all the pions from photons — a formidable task. On the other hand, with great experimental 
effort, it may be possible to count the neutral pions as well in the future. In that case, we will be able to distinguish 
DxC formation from other mechanisms of coherent pion productions. For example, one can count nt, the number of 
pions (both neutral and charged) in each rapidity window. For DxC formation, or any other mechanisms of coherent 
pion productions where the field is aligned with a random direction in the four-dimensional chiral space (ttx, ttj,, ttq, ct), 
the fiuctuation of rif is large. On the other hand, for mechanisms of coherent pion productions where the field is aligned 
with a random direction in the three-dimensional isospace (tTj, , tTj^ , ttq ) but without involving the a direction (i.e., the 
scenario discussed in Ref. it is straightforward to show that the fluctuation of nt is small. 

There are clearly mechanisms distinct from the DxC which can enhance pion number fluctuations beyond the 
most naive uncorrelated statistical estimate. For example, resonances produce more than one pion per decay and 
hence enhance fluctuations. Presumably these mechanisms lead to fluctuations which are characteristically small, 
with a variance ~ a/Af where a is of order unity (compare with Eq. (2.26), where the variance scales like {n)/Af, 
where (n) ^1). If the pion number fluctuations due to the DxCs discussed here are small, it will be difficult to 
distinguish DxC induced fiuctuations from fluctuations coming from other processes. It is therefore useful to have at 
least a qualitative estimate of background due to these competing processes. Some of these backgrounds have been 
calculated; for example, in Ref. for the kinematics of SPS collisions. We note, however, that such calculations are 
necessarily highly detail dependent. For our purposes, they are useful in that they may give an indication of the scale 
of competing processes. To see clear evidence of DxC production via fluctuations, one must be fortunate enough so 
that the characteristic scale of the DxC fluctuations (in some kinematical window) greatly exceeds fluctuations from 
competing mechanisms. Whether or not this is case, cannot be deduced reliably from theory given our present state 
of knowledge. Ultimately, this must be resolved via experiment. 

In summary, we have constructed signatures for DxC formation in heavy ion collisions which do not require counting 
the number of neutral pions. Instead, we suggest counting the number of charged pions produced, and a large 
fluctuation would be a signal of DxC formation. We believe these new signatures will be useful in searches for DxC 
at RHIC and LHC. 
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III. ENHANCING SIGNATURES FOR DISORIENTED CHIRAL CONDENSATES VIA CONDITIONAL 

PROBABILITIES 



In the previous section we discussed how charged pion fluctuations could be used as a possible signature of D^C 
formation. There remains, however, a serious potential problem in picking signal out of the noise if there is a large 
incoherent pion contribution in even the "best" kinematic windows. If it is possible, however, to reconstruct neutral 
pions reliably, one can enhance the signal-to-noise ratio by a judicious choice of data cuts. In this section, we discuss 
how this can be achieved. We study the distribution of R under the condition that the "first k pions" in our kinematical 
window are all neutral, where by the "first k pions" we mean the k pions with the lowest pT- In contrast to the usual 
result (R) = 1/3, the conditional expectation values are substantially shifted from 1/3, where the noises due to 
incoherently emitted pions arc peaked. As a result, one can greatly enhance the signal-to-noise ratio by making cuts 
to retain only those events with the "first k pions" are all neutral. 



A. Conditional Probabilities for a Single Domain 

Recall that we have, in Sec. II, calculated the probability distributions fo{go) and f±{g±) for a single domain of 
DxC. One can also study the probability distribution of the aforementioned variable R = nQ/{no + n±) = go/{go+g±) 
which represents the fraction of neutral pions among all pions emitted from the D^C. By reparametrizing distribution 



(2.15) in terms of i? = cos^ 6 and gt = go + g± = sin ip (the subscript t stands for total), one finds the following 



distributions. 



with the marginal distributions 



d^Pigu i?) - - . dgt dR, (3.1) 

t: \l R(l~ gt) 



dP{gt) = ft{gt) dgt = -J-^dgt, (3.2) 

dP{R) = MR)dR = -^dR. (3.3) 
2v R 

Note that while /o is drastically skewed towards the low end and f± is flat, ft is skewed towards the high end. 
This may sound counter-intuitive, but one must bear in mind that the emissions of neutral and charged pions a re n ot 
independent and there is no contradiction. Also note that fn is exactly as predicted in Ref. Jl|-^,|6|-p|j2^ (cf. Eq. (1.1)). 



We will rename this probability function as Fq{R), where the subscript "0" reminds us that the D^C is an isosinglet 
by construction. The distribution is plotted in Fig. 2a. It is obvious that the shape is qualitatively different from the 
Poisson-Gaussian distribution due to incoherent emissions. The expectation value of i? is 1/3, 

(i?)o= / RFo{R)dR^ 1/3, (3.4) 
Ja 

which has the simple interpretation that it is equally likely for the pion to be a ttq, 7r_|- or 7r_, and hence on average 
a third of the pions are neutral. 

This distribution is a consequence of the fact that we have assumed the D^C to be an isosinglet, a reasonable 
assumption on physical grounds. However, let's consider the distribution of R after the D^C emits a single ttq. The 
density matrix after the emission, which can be written as A ttq pdxC t^q, where A is a normalization constant, poxC 



is the density matrix defined in Eq. (2.14) and ttq annihilates a neutral pion. This new density matrix is not an 



isosinglet. It is straightforward to show that the probability distribution for this state is 



isin6lcos2 6'd6' ^R-R-^^^dR , , 

dP = ^ = (3.5 

/isin0cos2 0d^? J^lR-R-y^dR 



or equivalently. 



Fi{R) EE /(i?|lst pion is neutral) = RFo{R) j RFoiR)dR = ^R^/^. 



(3.6) 
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The distribution Fi{R) is plotted in Fig. 2b, which is drastically different from Fo{R). The distribution is skewed 
towards the high end, while Fo{R) is skewed towards the low end. Moreover, the expectation value of R is clearly 
pushed up: 

{R)i= [ RFi{R)dR = 3/5. (3.7) 
Jo 

So we have arrived at the intriguing conclusion that, if the "first pion" emitted from a isosinglet DxC is neutral, 60% 
of the pions subsequently emitted from the D^C are neutral, a huge enhancement from the original expectation of 
33%. 




This extraordinary statement certainly deserves more discussion. First, what is our criterion to decide which is the 
"first pion" ? The answer is simple: it can be any criterion. It does not matter as long as it is a priori equally likely 
to be a TTo, a 7r_|_ or a 7r_. The derivation just depends on our removing a neutral pion from the isosinglet DxC. It 
can be the first pion emitted in time, or the last one emitted in time, or even the 17tli emitted in time. The criterion 
can also be unrelated to the order of emission. For example, we can choose the "first pion" to be the one with the 
smallest polar angle. One can use any of these criteria to identify the "first pion" , and if it turns out to be neutral, 
then the R distribution of subsequent emissions is always given by Fi{R), provided we are in the large number limit. 
However, this is only true in this idealized scenario, when all the pions are coming from a simple D^C domain. In 
reality, some of the pions come from incoherent emission, and if the "first pion" turns out to be incoherently emitted, 
the expectation value of R of the remaining pions is still going to be 1/3, not 3/5. As a result, we want to choose our 
criterion in such a way that the "first pion" is likely to originate from the DxC and not from incoherent emissions. 
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Since D^C pions by hypothesis have low pr ~ 1/^, where L is the size of the domain, a natural choice is to use the 
pion with the lowest pt as our "first pion" . 

After clarifying the meaning of the term "first pion" , we move on to discuss the physical origin of the modification 
of the probability distribution of R. In a nutshell, we are seeing the physics of (iso)spin alignment due to Bose con- 
densation. To illustrate the point, let us first consider the following apparently unrelated Stern-Gerlach experiment. 
Consider a large number of massive spin-1 particles, which for concreteness will be called deuterons. Initially they 
are all polarized along a randomly chosen direction ft, which is a priori equally likely to be any direction in three 
dimensional space. In other words, 5 • n = for all the deuterons. Now let us pick one of these deuterons and pass it 
through a Stern-Gerlach spectrometer which measures Sz, the spin along the z-axis. What is the probability that the 
measurement gives Sz = 0? The answer is clearly 1/3, as the cases for Sz = +1, and —1 are equally likely. On the 
other hand, if the measurement on the first deuteron gives Sz = 0, what is the conditional probability for the next 
deuteron to pass through the Stern-Gerlach spectrometer also to be measured to have Sz = 0? The answer this time 
is no longer 1/3. The spins of all the deuterons are aligned along the same direction n, and that the first deuteron is 
measured to have Sz — suggests n is more probable to be more or less aligned along z than otherwise. As a result, 
the conditional probability is no longer 1/3, but can be easily shown to be 3/5, which is exactly the predicted value 



for {R)i in Eq. (3.7). The situation for a single domain of D^C is analogous, with isospin aligned pions instead of 
spin aligned deuterons. By construction, the pions in a D^C domain are isospin aligned, and by the same analysis, we 
have shown that the knowledge of the "first pion" being neutral can dramatically modify the conditional probability 
distribution of R. 

One can also consider the conditional probability distribution of R in the case that the "first pion" is charged. Note 
that 

Fo{R) = i(/(i?|lst pion is a 7r+) + /(i?|lst pion is a 7r_) + /(_R|lst pion is a ttq)), (3.8) 

and hence, since Fi{R) = /(i?|lst pion is a ttq), 

F{R) = /(i?|lst pion is charged) = \Fq{R) - \Fx{R) = |(1 - R)R-^/^. (3.9) 

The expectation value of i?, given that the "first pion" is charged, can be easily shown to be 1/5. As a consistency 
check, one can calculate (_R), the expectation value of R regardless of the species of the "first pion". Since the "first 
pion" is twice as likely to be charged as to be neutral, 

(i?)o = i(| + 2xi) = i, (3.10) 



agreeing with Eq. (3.4) 



Lastly, we will study the conditional probability distribution of R given that the k pions with the lowest px, which 
will be hereafter referred to as the "first k pions" , are all neutral. It is straightforward to show that in this case 



1; sin 9 cos^^ edO \R^ ■ R-^/^dR 
J isin0cos2'= ed9 ~ iRf^ ■ R-^/^dR 



dP^-Fi Ti -r^^ — (3-11) 



and 



Fk{R) = /(i?|lst k pions are aU neutral) = R^Fq{R) j R''Fo{R)dR = (fc + i)i?'=-i/2. (3.12) 

The distributions F2{R) and F3{R) are plotted in Fig. 2c and d, respectively. One can see that as k increases, the 
distribution is more and more skewed towards the high end. As a result, the expectation value of R increases with k. 

{R)k^ f RFk{R)dR^ {2k + l)/{2k + 3). (3.13) 

It is useful to define Q as the ratio of the number of 7r_|_ to the number of total pions emitted. By symmetry it is also 
the ratio of the number of 7r_ to the number of total pions emitted, and since R + 2Q = 1, 

(g)fc = l/(2fc + 3). (3.14) 

From the above analysis, the prescription to enhance the collective signal is quite clear. One should make successive 
cuts on the data sample on the condition that the k pions with the lowest pt are all neutral, and measure {R)k after 



each cut to see if it increases as predicted in Eq. (3.13). This result, however, depends on the assumption that we have 
only a single domain of DxC without any contamination due to incoherent pion emissions. Since this assumption is 
unrealistic for heavy ion collision experiments, the scenario we studied in this section is only an idealized situation. 
In the next section, we will discuss more realistic scenarios. 
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B. The Effects of Multidomain Formation and Incoherent Emissions 



The scenario considered in the last subsection is highly unrealistic in at least two ways. First, as discussed in the 
introduction, single domain D;\;C formation is highly unlikely. For a realistic treatment one must study D^C formation 
with more than one domain, each pointing in a different direction in the isospace. Moreover, we have neglected the 
effect of incoherently emitted pions, which have very important effects. If the neutral "first pion" is incoherently 
emitted, the R distribution of the remaining pions is described by fo{R), instead of fi{R) when the "first pion" comes 
from the D^C. In this section, we will incorporate these two effects and see how the predictions above are modified. 

We will study the expectation value of R, or equivalently the expectation of Q, for a situation described by the 
following parameters. The coherent fraction x is the fraction of pions which originate from D^C domains, so that 
when X = 1, all pions are coherently emitted, and when x — 0, all pions are incoherently emitted. We will consider the 
case where there are N domains, all containing an equal number of pion^, which will be assumed to be large. Each 
domain is described by an isosinglet density matrix, but the isospins of pions in different domains are uncorrelated. 
Now the question is: if the "first k pions" in this channel are all neutral, what are the expectation values of R and Q 
among the rest of the pions? 

The answer turns out to be the following expression: 

(i?> = i + 2A, (Q)==i-A. (3.15a) 

The shift A is given by 

A = xgP. - ^) ^ X(i - gP. (3.15b) 

where 

P, - p> (1 - pf-\ p = x/N. (3.15c) 

Each term in this formula has a simple interpretation: 

• The expectation value (R) is always 1/3 for the incoherently emitted pions. Only the pions coming from the 
domains are affected by isospin alignment; hence the outstanding factor of x- 

• Each coherently emitted pion comes from one of the domains, which will be called domain X. How many of the 
"first k pions" also come from domain X? The probability for each pion coming from domain X is p = x/N, and the 
probability that j of the "first k pions" coming from domain X is Pj — (^) p' (1 — pY~-' ■ 

• Given that j of the "first k pions" is coming from domain X, the conditional expectation value of Q decreases 
from 1/3 to l/(2j + 3), while the conditional expectation value of R increases by twice the above quantity. 

In passing, we note that A can also be expressed as an integral or the hypergeometric function 2^1: 

,l^cos^^^3e 

A3 sin^e Jo cos2'=+4^y' 



X 



I- x/N 

(3.16) 

In Fig. 3, we have made contour plots of A = 1/30 (such that {R) — 0.4 and (Q) — 0.3), for k — 1, ... ,5 in the 
(x, 1/.^) parameter space. The horizontal axis is the coherent fraction x while the vertical axis is 1/A^ where N is 
the number of domains. Both x and 1/A^ range from to 1. Thus, for example, with 3 domains and x = 0.6, in order 
to have A > 1/30 we must have fc > 3. 



Equations (3.15) illustrate the main results of this section. One can see that, without any DxC formation, x — 



(corresponding to the left edge of Fig. 3), A vanishes, and {R) = (Q) = 1/3 as expected. The bottom edge of the plot 



^This assumption of all domains having the same number of pions is unrealistic but is made for illustrative purposes. The 
effects of unequal domain sizes will be briefly discussed below. 
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corresponds to iV ^ oo and also gives A = for any finite value of k. The shift A is largest for a single domain of 
D^C without any noise due to incoherently emitted pions, i .e., w hen x = TV = 1 (the top right corner of the contour 
plots), giving A = 1/3 — l/{2k + 3) and reproducing Eqs. ( 3.13| ) and ( 3.14 ). For fixed values of (x, N), A increases 
with fc, accounting for the spreading of the parameter space with A > 0.4 as k increases from 1 to 5 in Fig. 3. 




FIG. 3. Contour plots of A = 1/30 (such that (R) = 0.4 and (Q) = 0.3) for different values of k. The horizontal axis is the 
coherent fraction x, while the vertical axis is 1/A'^ where A'^ is the number of domains. Both x and 1/A'^ range from to 1. The 
curves, from top right to bottom left are for fe = 1, 2, 3 4 and 5, respectively. The shift A is larger than 0.4 above the curves 
and smaller than 0.4 below the curves. 



One expects that when the number of DxC domains is large {N ^ 1) or when most of the pions are incoherently 
emitted (x ^ l)i it will be difficult to observe clear signals of DxC formation. However, in such situations x/^ is 
small and A is dominated by the j — 1 term (the j = term always identically vanishes) and 



2x^k 
15N 



Oi 



(3.17) 



Thus a large k may make up for a small coherent fraction Xi or a large number of domains TV, and enhance A, 
whic h des cribes the shift of (R) and (Q) from 1/3, to an experimentally measurable magnitude. From the form of 
Eq. (3.17), one expects this shift to be substantial whenever k ^ N/x^- However, even for a value of k as small as 
TV/4x^, A = 1/30 + Oix^/N"^), which translates to (i?) = 0.4 and (Q) = 0.3 — a substantial deviation from the 
incoherent case. This suggests one should make successive cuts for events where the k pions with lowest px are all 
neutral, and study (i?) after each cut. An increase of (i?) with k would suggest that DxC domains are formed. 

Equation (3.17) appears to suggest that one can increase A to an arbitrarily large magnitude by choosing a 
sufhciently large value of k. Of course this is not true. Equation ( 3.17 ) is obtained as the leading term in a x/^ 
expansion, but when k oo, this expansion breaks down as terms of higher order in x/^ ^^'^ enhanced by factors of 
(^.). We can easily see that 



^X, (i?) ^ i + fx, 



(Q) ^^(1 — x); fc^oo with X and N fixed. 



(3.18) 



In other words, our signal enhancement scheme is fundamentally limited by the amount of noise due to incoherently 
emitted pions. When x is small, most of the pions are incoherently emitted, and for them, (R) is always around 1/3 
regardless of what cuts one makes. On the other hand, the large k limit of (R) does not depend on N , th e number 
of DxC domains. Recall that we have several distinctive signatures, like the R distribution in Eq. (1.1) and the 
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conditional expectation values for R described in the previous section, for a single domain of D^C, where all the pions 
in the D^C are isospin aligned. With multi-domain formation, where the pions in different domains may point to 
different directions in isospace, the effect of isospin alignment is greatly washed out. However, given that the "first k 
pions" are all neutral with fc ^ A^, it is probabilistically extremely likely that each of the N domains is the origin of 
some of these "first k pions". As a result, each of these N domains are well- aligned along the ttq direction, and hence 
also well-aligned with each other. As fc — > oo, the N domains look more and more like a single big domain in the ttq 
direction, which is the case where the signature is the most dramatic. In short, with large k, our cuts are picking out 
the events where the signals are the strongest, and hence resulting in a large signal-to-noise ratio. 

In above we have assumed the sizes of all N domains are identical for illustrative purposes. A more realistic 
treatment would have N domains, each with different sizes pi, 1 < i < iV, such that 'Y^^Pi = X- (The size of a domain 



is defined to be the fraction of pions which originate from this particular domain.) Then Eqs. (3.15) are generalized 
to 



^ 1 
3^ 



N 
1=1 



fc 

E 



k-j 



1 



2i + 3 



(3.19) 



In the weak signal limit, i.e., when all the piS are small, Eq. (3.17) gets modified to 



2fc 
15 



N 



J2p^ 



i=l 



^3 



2k 
15 



P + Oi 



7V2. 



(3.20) 



where p = '^^pf has the following nice interpretation: p is the average over all pions (both coherently and incoherently 
emitted) of the sizes of the originating domains, which is pi fo r a p ion from domain i and zero for an incoherently 
emitted pion. For N domains of equal sizes, p — 'x^ /N and Eq. ( 3.17 ) is recovered. Again we s ee that A grows linearly 
with k in the weak signal limit. As fc ^ oo the shift A is again limited by the bound (3.18), which applies also for 
the cases of unequal domain sizes. 



To recapitulate, we suggest the following experimental procedures: 

• Count the number of neutral and charged pions event by event from heavy ion collision experiments and measure 
their individual transverse momenta and rapidities. 

• Apply a low pt cut to suppress the noise due to uncorrelated pion emission. 

• Bin the events in different rapidity windows. 

• In each rapidity window, calculate the expectation value {R). 

• Make a cut to retain only events where the pion with the lowest pT is neutral. 

• Calculate, in each rapidity window, the expectation value {R) for all remaining pions in all events which survive 
the cut. 

• Make another cut on the surviving events to retain only those where the pion with the second lowest pT is also 
neutral. 

• Again, calculate in each rapidity window the expectation value {R) for all remaining pions in all events which 
survive the cuts. 

• Repeat the above prescription of making successive cuts to retain only events in which the pion with the next 
lowest Pt is also neutral, and calculate (i?) for each rapidity window after each cut. If we find (i?) deviates from 1/3 
then we are seeing signatures from D^Cs. 

Note that this prescription requires reconstructions of pt's of individual pions, both charged and neutral. We have 
also presumed that the coherent fraction x and the number of domains formed N are roughly the same for each event. 
(More specifically, the probability distributions for x and A'' are narrow peaked.) 

By applying these successive cuts, we are retaining the events with D^C formation and most of the pions are 
well-aligned along the ttq direction. What is being cut are the events with D^C formation but most of the pions are 
well-aligned along the t:^ or tTj, directions, and the events where there are incoherent pions with very low pT, which is 
the main source of noise to our signal. As a result, these successive cuts are substantially improving the signal-to- noise 
ratio, making it easier to observe D^C formation. On the other hand, just like any other cuts on data to suppress the 
noises, we are giving up on statistics. Moreover, for large k we are cutting on rare events so the loss in statistics can 
be severe. For the cases where the signal is weak (small coherent fraction x <C 1 or large number of domains 1) 
on each cut we are losing about two-thirds of the events. 

In conclusion, we have devised new cuts to enhance the signal in searches for D^C. These cuts retain only events 
where the k pions with lowest px are all neutral. We have shown that, after these cuts, the fraction of neutral pions 
within the remaining sample is substantially larger if D^Cs are formed in the heavy ion collision. 
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IV. SUMMARY 



In this paper, we have discussed two methods for detection of D^C in heavy ion coUisions. Both methods make 
use of the fact that all the pions in a D^C are aligned in the 0(4) chiral space, and hence also in the 0(3) isospace. 
This leads to a larger fluctuation in the number of charged pions produced than the case where all the pions are 
independently emitted. As a result, one can look for signatures for D^C formation from charged pion data only. The 
alignment in isospace also suggests that, if the "first k pions" are all neutral, the remaining pions are also much more 
likely to be neutral than otherwise. By making cuts accordingly, one can greatly enhance the signal-to-noise ratio. 
We expect these methods will be useful in future searches for signatures of D^C formation in RHIC and LHC. 

Support of this research by the U.S. Department of Energy under grant DE-FG02-93ER-40762 is gratefully ac- 
knowledged. 
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